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ABSTRACT 


The range of fluid dynamics applications in the lung 
is vast, encompassing many interesting problems rarely 
touched in any other single engineering area. A survey 
of recent work in this field was conducted. Additionally, 
a numerical solution for gas mixing in a simple model of the 
lung was developed indicating its application to a wide 
range of similar problems. Finally, new parameters were 
experimentally investigated in flow оро a thin-walled 
collapsible tube. Results indicated the effects of .the 
parameters to be significant in the flow of blood in the 
pulmonary capillaries as controlled by the surrounding 


aveolar pressure. 
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I. SURVEY OF FLUID DYNAMICS IN THE LUNG 


A. INTRODUCTION 

The study of fluid dynamics in the lung is a wide field, 
encompassing a large range of problems. The vastness of the 
field is easily understood when the complexity of the organi- 
zation of the respiratory tract and the cardiovascular system 
is considered. Many of the problems which arise in these 
areas are unique even when the fluid dynamics of the problems 
wec ons idercedmn 1solation. Often, these problems require 
work in areas of fluid dynamics that have never been 
Merenpted in purely engineering applications. Frequently, 
new researches in the field center on very basic fluid flow 
problems, and the results are as interesting as any that are 
mepested Dy engineerine applications. Figure l indicates 
the main areas of study within the cardiovascular-respiratory 


system. 


B. FLUID FLOW IN LARGE AIRWAYS AND ARTERIES 

In both the se airways of the lung and in the larger 
arteries the Reynolds number may vary from several thousand 
to one or two hundred. This range of flows has not been 
ellosely studied, particularly the large Reynolds number end. 
In most engineering applications large Reynolds number 
studies have been limited to fully developed turbulent flow. 
In the lung, however, high Reynolds number flows are only 
sporadically or locally turbulent, seldom becoming fully 


developed. 





Inn studies Буена аша Ostrach [Ref. 1] for large 
Reynolds number flow regimes including the effect of 
tic cosity, pulsatile flow, and vessel elasticity, the 
developing flow was found to have a strong boundary layer 
character, resulting in high wall shear stresses near the 
vessel entrances. The development length was found to be 
large and dependent on the Reynolds number and unsteadiness 
parameter, while wall distensibility was not significant. 
Thus flow development characteristics are ЕЕЕ сапе in the 
large airways and arteries. 

Additionally, the large vessels may be subject to 
curvature, inducing large secondary flows due to centrifugal 
forces. These forces, acting as a result of curvature of 
the primary flow streamlines, produce secondary flow 
velocities at right angles to the primary streamlines. 

Study by McConologue and Srivastava [Ref. 2] has indicated 
that the resulting motion of a central fluid is toward the 
outside of the bend, and a return flow toward the inside of 
the bend near the wall. The motion toward the outside wall 
15 nearly uniform while the return flow, retarded by viscous 
stress, is limited to a thin boundary layer. A result of 
this fact was the increase in critical Reynolds number 

based on pipe diameter for transition in curved pipes 
compared with straight pipes. In addition, the secondary 
flow promotes a much more rapid lateral mixing than would 


be possible for normal laminar flow. 





Work has been undertaken to study the combined effect 
БІЗ the last two topics апе forthe fact that the first entry 
region to both the cardiovascular system and respiratory 
tract is a curved tube: the great arch of the Aorta, and 


the mouth and trachea. 


C. COMPLEX BRANCHING 

Complex brane hane oi tene passageswexists ІП БӘСІ (пе 
respiratory tract and the cardiovascular system. As the 
flow branches early in each system it undergoes a large 
distortion at each branch, and the flow may not recover 
before reaching the next bifurcation. In the smaller air- 
ways and capillaries where the Reynolds number is too small 
for this effect, the flow becomes difficult to analyze 
due to the large number of tubes in parallel and the lack 
of information on the distribution of tubes by length and 
diameter. 

Pedley, Schroter, and Sudlow [Ref. 3] found that in 
the higher Reynolds number airways secondary flow was 
exhibited downstream of each bifurcation, with peak axial 
incl ties Just outside the boundary layer оп the inner 
wall of the junction. These secondary motions decayed 
slowly downstream. 

Additional study by Caro, Fitz-Gerald, and Schroter 
[Ref. 4] has indicated that the rate of shear in a branching 
system was dent on two factors: a flow factor depend- 


ing on the total flow in a particular tube, and a geometrical 





factor. The geometrical factor was found to decrease 
downstream. It was also greater on the outside of a bend 
than on the inside. At a bifurcation the factor was largest 
on the central flow divider, the outer bends of the two 
temeanter branches becoming inside bends, effectively. The 
flow factor was also found to decrease with distance down- 
stream. The general result was that the rate of shear was 
least nearer the heart or mouth and, within these, at the 
meee Where the geometrical factor was smallest. 

In steady flow a branching system affects pressure 
Ggaeps at each series of bifurcations. Here the effect of 
resistance differs for the arterial problem and the 
bronchial problem, as discussed by Lighthill [Ref. 5]. In 
the bloodstream the greatest pressure drops occur in the 
smaller vessels of diameter less than 0.1 mm., whereas in 
the airways the maximum pressure drops occur early in the 
branching, soon after the fifth generation. The airways 
maintain a much larger daughter-to-mother branch ratio 
CHOSET #6 the system entrance. This coupled with the higher 
Reynolds number in those regions would tend to cause separa- 
tion, a possible explanation of the phenomenon. The arterial 
resistance problem is said to be a low Reynolds number 
problem while bronchial resistance is considered a high 


Reynolds number problen. 


D. BLOOD MICROCIRCULATION 
Flow in the small capillaries or the lung differs 


greatly from that usually found in engineering fluid 





dynamics. The feature which distinguishes this type of 
Blowes that шем иташшйеттта effects are itotally negli- 
gible. The Reynolds number is so low that the flow pattern 
is determined solely by the need for pressure gradients and 
viscous forces to be in exact equilibrium. Under these 
conditions, Bernoulli's equation is not applicable and there 
is essentially no difference between static and dynamic 
pressures. Centrifugal forces are also negligible: fluid 
can negotiate sharp bends with no secondary flow effects. 
Ші сепета!; fluid motion is insensitive to vessel geometry 
and there is little tendency for separation. This insensi- 
tivity to geometry was further demonstrated by Fitz-Gerald 
ПЕ. 6] with respect to entrance length effects. fhe 
difference between Poiseuille motion of a fluid and fluid 
entering a tube at a large angle, with the same total flow, 
was reduced by an order of magnitude in an axial distance 
of half the tube radius. 

Pro pasao canes рае сове ene large vessels 
governed mainly By a balance between distensibility and 
inertia, 15 governed in the smaller vessels entirely by a 
balance between distensibility and resistivity. A feature 
of the blood circulation in most parts of the body is the 
vasomotor control of peripheral perfusion, the ability of 
smooth muscle to alter the diameter of small vessels in 
order to vary the flow rate through the vessels by an order 
of magnitude. Contrary to the systemic microcirculation 


system, microcirculation in the lung exhibits almost по 


10 





vasomotor control in the pre-capillary vessels or elsewhere. 
As noted by Lighthill [Ref. 5], most of the resistance to 
flow appears to be in the capillaries themselves. 

Perfusion in the capillaries can vary greatly with 
location in the lung, an effect based predominantly on 
Мат су as presented by Glaister [Ref. 7]. A typical 
pulmonary capillary is easily distorted by pressure dif- 
ferences between the capillary and alveolus. In effect, 
the capillary is a thin-walled collapsible tube in which 
flow is governed by the alveolar pressure Р, and the 
upstream and downstream capillary pressures, Р, апа P» 
respectively. The following information was obtained and 
Aai fied by Conrad [Ref. 8], Holt [Ref. 9], and many 
others. The effect of gravity was to divide the lung into 
иисе basic zones. In zone l, at the top of the lung, Р, 
was greater than Р, or P^, and there was no flow. In zone 2, 


the middle of the lung, Pj was less than Da which was less 


2 
than P1» and flow was governed by the pressure difference 
Pı"P,- In zone 3, low in the lung; 9% was less than Р, or 
РЬ, апа flow was unobstructed and proportional to the 
pressure difference Ру-Р„. An experimental investigation of 
flow through a thin-walled collapsible tube is included in 


this paper in Section III. 


СЕК FLUID EROPERTTES 
A problem relatively unique to physiological fluid 
flow lites in the unusual fluid properties of blood, and to 


a lesser extent of air. Blood is a suspension (40-505 by 


11 





volume) of small deformable bodies, mainly red cells, in a 
transparent fluid, the blopd plasma. Although the flow 
properties of plasma are similar to those described by 
Newton's viscosity law, the concentration of red cells at 
low strain rates yields substantially increasingly effective 
viscosities. However, Lighthill [Ref. 5] indicates that 
data suggest that at the higher strain rates in the larger 
vessels the errors are not large in treating the blood as 

a Newtonian fluid with constant viscosity. Studies of 
particulate suspensions by Mason and Goldsmith [Ref. 10], 
simulating blood over a range of Reynolds numbers extending 
from creeping viscous flow to that where inertial effects 
play an important role, indicated that at low Reynolds 
numbers for a single particle, particle deformation and 
interaction with the wall forced migration to the axis of 
the tube. At higher Reynolds numbers, particles also 
migrated laterally, but to an eccentric radial position. 

In blood-similar suspensions at low Reynolds numbers, 
blunting of the velocity profiles was evident, the dispersions 
were РЕ ЕЭ and cells were deformed significantly. 
At higher Reynolds numbers in steady flow, and in oscillatory 
and pulsatile flow, cells tended to migrate inward from tne 
tube wall. This resulted in a two-phase flow of a central, 
more viscous core surrounded by a thin, particle-depleted 
and less viscous zone, with an attendant decrease in the 


energy dissipated in the tube. 


12 





AppTtieatioms#eor ndreodynamac- lubrreatıionrtheory!" Бу 
ВИЕ Ш [Ref. I1] to blood flow in the narrowest- capil- 
laries have corresponded to the observed facts that the 
stress necessary to deform a capillary wall was nearly two 
orders of magnitude larger than that required to deform a 
red cell. Thus it is that the cells ШТ very substan- 
tially in order to pass through the smallest capillaries. 
Lubrication theory further suggests that the stress, which 
the wall of a very narrow capillary passage must transmit 
a red cell sliding past it so that the cell will deform 
sufficiently, must be transmitted by viscous action in a 
thin lubricating layer of plasma; that the plasma film 
thickness must adjust passively to the value required; and 
that sliding may cease when the plasma film thickness 
becomes too small. Theoretical observations also indicated 
that the film thickness should vary with the square root of 
the cell velocity. This is an interesting non-linear effect. 

Air is also a suspension of dust and other fine 
pieles, and significant study has been devoted to the 
behavior of airborne particles in the lung. The deviation 
of particle motion from gas motion is caused by Brownian 
motion, gravity, and inertial forces. Altshuler (Ref. 12] 
Stated that for particles of less than 0.2 ит the dominant 
force is Brownian motion. For particles of greater трап 1 уш 
the dominant effects are settling and inertial displacement. 

Owen [Ref. 13] observed that for inspired flow rates of 


40 L./min. and greater the boundary layer on the tracheal 


13 





wall was ange Ее іт ессе, 1f- the air 

entering the trachea separates at the vocal chords, almost 
fully developed turbulent pipe flow may exist and decay 
slowly. This would suggest that many particles are deposited 
on the tracheal walls; however, the normal rate of evapora- 
tion from the tracheal wall causes a radial motion sufficient 
to impede the deposition of all but the finest or coarsest 


punticles. 


G. GAS MIXING IN THE RESPIRATORY TRACT 

дето са breathing tregkiencies the main mechanical 
impedance to gas flow within the lungs stems from cu elastic 
distensibility of the air spaces terminating the airways 
(alveoli), the flow resistance of the air passages being 
ener emely small. Due to the weight of the lungs, there is 
Weradient of increasing distension from bottom to top. 
Mead [Ref, 14] found that the lower, less distended regions 
шетте more compliant than the upper, more distended ones, and 
that there was a gradient of increasing gas flow from top to 
bottom. As a result, the ОСБ оспе lungs received 
ОШГУ half the bulk flow of gas per gram of lung tissue 
received by the lower regions, 

Davidson and Fitz-Gerald [Ref. 15] discovered that for 
various geometrical models for small airways in the lung, 
the bulk flow was significant for gas transport in the final 
"711127 thoa lung: 

In single-breath nitrogen washout experiments, as 


discussed by Farhi [Ref. 16] and Pedley [Ref. 17], in which 


14 





the subject inspired a single breath of pure oxygen, then 


рге Iene N Concentration was mneasuredeat the mouth 


2 
daning expiration. eine results indicate that initially 
БИЕГЕ” 15 ПО N, because the expired gas is the pure 0, most 


decently inspired, which never made contact with the N, 


rich air already present in the lung before 0_ inspiration 


2 


Began. Following this region, the N5 concentration rapidly 


2 
imereased to a nearly constant value on the "alveolar 
peetecau."' The alveolar plateau represented air from deep 


апе luneszwhere”the N- concentration was initially that 


2 
of ordinary air. 
Gas mixing in a simple model of the lung was solved 


Шинг-г1Са11у in Section II of this paper. 


G. GAS EXCHANGE AT THE ALVEOLAR INTERFACE 

The interchange of materials (05, C05) Бу йе Circulatory 
system depends on two processes: (1) circulation of blood 
through the capillary beds and (2) transport of various 
components of the blood through the capillary walls. Net 
transport across capillary walls is driven by two basic 
ТЕБ: (1) differences in hydrostatic or effective 
osmotic pressure of the fluids on either side (ultrafiltra- 
tion transport) and (2) differences in concentrations of 
permeating solutes on either side (diffusion transport). 

Ren tii Rer ШЕ слона си сола оваа 1 Стат 1 traitdon 
transport was brought about by changes in the ratio of the 


щаба ватсещо вее невеста ща е Of 2 capillary bed to the 
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total resistar e ol the Бей, апа Б спапсев in the number 
of open capillaries. Control of diffusion transport was 
broucht abent by changes im the number. of open capillaries, 
and by changes in local distribution of capillary blood 
flow relative to capillary surface and permeability. 

‘ As blood flows along the pulmonary capillaries it 
meeeiyes oxygen from the alveoli, and delivers carbon 
Тас то Lich Пато concentration gradients for the gases 
wevelop along the capillaries, Increasing for oxygen and 
decreasing for carbon dioxide. Study of digital computer 
models of pulmonary gas exchange by Milhorn and Pulley 

E 19] revealed that the partial pressure of blood reached 
equilibrium with alveolar oxygen at about two-fifths of the 
Geral distance along the capillary pathway. Also the partial 


pressure of the CO, reached equilibrium with the blood about 


2 
halfway along the capillary pathway. Additionally, the 
Buatral pressure gradient of oxygen had a concave shape due 
HChe Convex shape Of 1tS dissociation curve. The concavity 
dd not exist in er Pencon entrato pradient of oxygen 
or in each of: the curves for carbon dıoxıde, 

Ine eitect of pulsatile flow of blood on oxygenation of 
the blood in pulmonary capillaries for an ideal lung was 
investigated by Graeser, Kim, and Crandall [Ref. 20] for 
several combinations of exercise and atmospheric oxygen 
Content. The pulsatile flow was found to give substantially 


the same end capillary alveolar-arterial tension gradient 


as that found for a constant velocity model. When breathing 


16 





air, the pulsatıle flow effect yielded a 10-15% greater 
dir fusing capacity thon the constanti velocity model, bùt 


no effect was observed in low oxygen cases. 
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II. NUMERICAL SOLUTION FOR GAS MIXING 


IN A SIMPLE MODEL OF THE LUNG 


A. INTRODUCTION 

The diffusion of gas throughout the lung is a complex 
problem dependent on both convective transfer and molecular 
ч ти<1оп. A typical physiological experiment in this field 
15 the single breath nitrogen washout experiment described 
by Farhi [Ref. 16] and Pedley [Ref. 17]. The subject, 
Sear ting from an arbitrary lung volume, inspired a single 
breath of pure oxygen. The subject then exhaled and the 
№, concentration at the mouth was measured during expira- 
tion. The resulting graph of N, concentration versus time 
was typically similar to that of Figure 2, in which high- 
frequency oscillations due to the cardiac cycle have been 
averaged out. In region I there was essentially no N, due 


2 


to the fact that the expired gas was the pure 0, most 


recently inspired, which never reached the N, а О а 
already in the lung upon inspiration. Region II was a 
transition region where the N, concentration increased 
rapidly with time, reaching a nearly constant level called 
the "alveolar plateau" in region III. The alveolar 
plateau represented air from deep in the lungs where the 
initial N, concentration was that of ordinary air. If the 
total volume of 0, inspired am (related to time by the 


mass flow rate) was less than the anatomical dead space (Vp) 


(the volume of the airways as opposed to the alveoli), then 
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the volume Vo at which the plateau was reached was 
approximately equal to ee: Ши NES was greater than Ур» 
ү 


g was approximately equal to V although the exact relation- 


р? 
ship was a function of the gas mixing and diffusion of the 
inspired gas with the air previously present. 

A significant effect in region III was that the N, 
concentration was not constant, but exhibited a marked slope. 
There are two possible explanations for this phenomenon. 

The first is based on "regional inhomogeneity," that is, 


the concentration of N, in the alveoli near the top of the 


2 
mung is less than at the bottom because the upper regions 
we initially more expanded than the lower due to the 
difference in pleural pressure. Thus, less oxygen may 
enter the lower regions during inspiration. This also 
requires that the air from the upper regions is expired 
Hst, often called “sequential emptying." There is evidence 
supporting regional inhomogeneity; however, sequential 
emptying is quite controversial. The second explanation is 
that of "stratified оше ое Сша iS athere is a 
concentration gradient between the inspired O, and the N, 
rich gas already present, and there is no a priori reason to 
assume this gradient will not exist in the alveoli. 

The simple model en follows was developed to study 


the gas mixing problem, and the N, washout case in particular. 


2 


19 





Ё 11 Т 2! Р? 

The basic development of this model was suggested 
by Pedley [Ref. 17]. 

It has been noted that the lung is a very complicated 
network of branched tubes, іп which the exact velocity 
distribution remains unknown. Reynolds numbers may vary 
from several thousand to less than one. Thus, solution of 
Pmetull dittusron equation, for which the fluid velocity 
weld is required, is impossible, and considerable simpli- 
fication is necessary. 

The lung was assumed to be a symmetric, dichotomously 
branching system, in which all pathways have the same 
length. The lung was replaced with a one-dimensional 
system whose cross-sectional area a(x) increased rapidly 
from the trachea to alveoli as shown in Figure 3 and 
consistent with Weibel's [Ref. 21] data. The area was 
considered constant from the mouth to the first branching 
of the airways, where it increased greatly. The cross- 

2 


sectional area of the trachea, a is tapproximately 2.5 Cm 


0? 
The area of the final generation of airways is of the order 
of то см“. The position where the area began to increase 
was taken to be x = 0. The distance from the mouth to 

Х - 0 is about 38 cm, while the remaining distance to the 
final generation 1S about 6.6 cm. a(x) was assumed to be 


a continuous smooth function independent of time. А11 


volume change was considered to take place in the alveoli. 
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The N, concentratıon was considered to be uniform across 
JN NUSS section. “ГС іле average concentration c(x,t) 
was a function only of distance x down the lung and of time t. 
The velocity u(x,t) was also assumed to be uniform across 
a cross section. These assumptions were not as drastic as 
pa first appear: The trumpet-shaped tube in Figure 5 15 
actually a large number of small tubes nearly in parallel, 
e pecially at thè alveolar end. The main effect of radial 
са ст слез velocity profiles, and turbulenee isto 
лепсе the relationship of concentration to the longi- 
tudinal velocity. Compensation for this omission was 
imenoduced іп the diffusivity. 

Previous theories by Cuming, Crank, Horsfield, and 
Parker [Ref. 22] and La Force and Lewis [Ref. 23] have used 
similar models, but greatly oversimplified the actual gas 
mixing process. They ignored the fact that mixing takes 
place everywhere simultaneously with convection, especially 
in the upper airways. They assumed a period of convection, 
in which the ЕАН ретуееп риге 0, and mixed air remained 
flat, ОКОТО E TO е сув ол ома T0llowed by another 
PDerrod of convection. Secondly, they assumed that the 
dii fusion was governed at all times by the molecular diffu- 
sivity Dol’ Ша reality оо laminar and turbulent flow, 
ene existence of a velocity profile or turbulent eddy ensures 
that longitudinal mixing is much greater than can be described 


by püre molecular diffusivity де Ток [Refs. 24 and 25] found 


that after a certain time the average concentration could be 
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calculated from a knowledge of the average longitudinal 

velocity, on the basis of an effective longitudinal 

an fusivity mach greater than Dol’ 
Flow in the trachea was assumed to be turbulent and 

an effective longitudinal diffusivity (Dg) was calculated 

using Taylor's relation. For a tracheal Reynolds number of 


5000, D, was approximately 2060 cn^/sec. In the terminal 


0 
regions of the lung, the mixing is governed almost entirely 
by the molecular diffusivity with a value of about 0.2 сп2/ 
sec. Іп between these values, the flow undergoes large 
changes, about which little is known specifically. .D(x) 


was considered to vary smoothly between D, and D from 


0 ol 
и И со the alveolar ena, Therefore, the important 
feature of greatest mixing at the highest velocities was 


introduced. 


C. ANALYSIS 
The governing equation for diffusion in the one- 


dimensional model-was: 


әс дс 9 оС D 9a дс 
iwc a CD И a Oe Ox 


where the last term accounts for the change in cross 
sectional area. 


Prom Continui ty: 


и „а 


uc t) = aT 
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where ug was tes locı in the trachea, and ag was the 
cross-sectional area of the trachea. Inspiratory and 
expiratory flow rates were considered uniform and to be equal 
and opposite. 

The boundary conditions on c(x,t) were as follows. 
Both с and дс/дх were continuous at x » 0. At the alveolar 
end, for N, there was no flux across the end, 9c/ax = 0. 
At the mouth, c = 0 for inspiration (pure 0,) апд дс/дх + 0 
for expiration (diffusion negligible compared with 
convection). | 

The initial condition for inspiration was uniform 
Eouccntratiron everywhere. The initial condition for expira- 
пощи наз the distribution of concentration at the end of 
inspiration. 

In the governing equation for diffusion, the concentration 
was non-dimensionalized with respect to its initial value, 


and all other quantities with respect to their values in 


эв 


the trachea (C = D = ЭЭ X E ES S Ше Би? 
: D Кы? qup & 


T 0 0 


The following equation was obtained: 


22 


9Х 


e) 


ae 9 Mo Mee 
эт AN + [+N - 57 (AD)] zx = AD 


г 


where the minus sign denotes inspiration, the plus sign for 
expiration, andN = 40/00» the effective Peclet number in 


tne trachea. 
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The new boundary and initial conditions were: 
ze! 
OM: о раа OS 


QS TE E 


34 (1-6,1) «07 


for expiration: 


3C _ ac z 
5% С-5,7) < << (1-6,7) = 0 


where 6 = a.) апат иес отав о тсасбопаБ е 
Peclet numbers, the values of A and D were close to their 
actual values (8 2 0.02), at X = 1-8. 


The particular forms chosen for A(X) and D(X) were: 


SEN S 0 A= Del 
ЖОО NOS) А - (1-Х)77 

= mol 
17: D-p ,+—3 


the function of A = (12547: was consistent with Weibel's 


data. 


D. SOLUTION 

Solution of the non-dimensional equation was performed 
numerically with the use of finite difference approximations: 
ЗС/9Т by forward difference, and 9С/9Х, 3^c/ax* by central 
difference. The resulting equation was solved on the digital 
ОРЕСТ ИБА е перете еу оша time aind distance "Steps Со 


maintain stability and accuracy. The "concentration at any 
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tine T was calculated from data for the previous time step. 
Solutions were generated for four values of the effective 
Peclet number in the trachea. Inspiration time before 


expiration was also varied. 


БКО RESULTS ANB DISCUSSION 

The simplified model of the lung produced results 
Ear to those described by Farhi [Ref. 16] and Pedley 
[Ref. 17], with a minor variation. In addition, several 
interesting results were observed. 

The effective Peclet number in the trachea was found 
to have an effect that varied with distance down the model, 
as shown in Figure 4 and by data in Appendix A. In the 
trachea, an increase in the Peclet number caused a decrease 
in the N, Concentration at a particular position. This 
ertect was reversed for all X greater than zero, in the 
actual branching. That is, an increase in Peclet number 
increased the N, concentration at a particular position. 
In addition, this -effect tended to be felt just upstream of 
X = 0. The general result was that as the Peclet number 
decreased, the further down the lung the diffusion progressed. 
This indicates that a reduced Peclet number corresponded to 
ШОЧ туе increase In the diffusivity with respect to the 
velocity. 

Diezerrecteosssene duration of inspiration was also 
dependent on distance down the model. In the trachea, the 
concentration rapidly decreased with time, while in the 


lung the process was much slower, as shown in Figure 5. 
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N, concentration in the model during expiration also 
varied considerably from the trachea to the lung. As 
epi ration time increased, the N, concentration slowly 
increased linearly along the trachea, then rapidly increased 
with distance to the alveoli (Figures 6, 7, 8). Increasing 


expiration time also produced an increase in the N, concen- 


р 
tration at a particular position, with the exception of the 
last few airways before the alveoli. 

lie N, washout curves obtained were similar to those 
pected with the exception that the sharper changes in 
ш рс were not detectable between regions II and III, as 
шиш in Figures 9, 10, and 11. In fact, for the longest 
inspiration (T = 96), the curve for Peclet number of 12 
was essentially linear. For this reason it was very 
difficult to examine the relationships between Vg» Vp» 


and Mr as discussed in the introduction. Үр was 55 


(non-dimensional), and for Figure 9, V, = 12, and the 


0 
relationship appeared LOM MOG Е ОЕ а ое ОУ, ШЕ = Vo = 48, 
but it was difficult to determine exactly where the alveolar 
plateau began. Бод ариг? Vo = 96 and the alveolar 
plateau should have begun at Vp = 55. It obviously did not 
for either case. 

The effect of Peclet number was to decrease the 


concentration at a particular time with an increase іп 


Peclet number for the washout curves. 
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F. CONCLUSIONS 

The results of the numerical solution for the 
theoretical simple model of gas mixing in the lung justified 
the following conclusions: 

1. The model effectively simulated gas mixing in the 
lung with simultaneous diffusion and convection БІРЕСЕ, 

2. The model is very flexible and many parameters may 
be easily varied to more closely approximate the experimental 
data for the washout curves. 

БР "с келен ity Of the model extends to many other 
Similar experiments in the field. The ease with which 
initial conditions, boundary conditions, and the various 
parameters may be varied and the rapidity of computer 
solution provide a fast, accurate tool for future research 


in the field. 
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БЕП 20,455 ай ОО ИО СЕВО 


IN A THIN-WALLED COLLAPSIBLE TUBE 


A. INTRODUCTION 

nvas noted earlier, іп the suryey Section of this 
paper, that flow in the small capillaries of the lung is 
analogous to flow in a thin-walled collapsible tube. Flow 
ШІ (115 type Of tube 15 dependent on the surrounding pressure 
as well as the upstream and downstream pressures. 

The force of gravity exerted on the lung tends to 
divide it into three basic zones. In zone I at the top of 


the lung, the alveolar pressure P, is greater than the up- 


А 


Stream static pressure P in zone II, near the middle of 


17 
1s greater than P 


pae lung, P which is in turn greater 


1 А 
than the downstream static pressure P. ША” theo est 
n 15) lessmenan both Р, апа 


P». The pressure-flow relationships in these three zones 


por the lunp, zone III; P 


have been experimentally studied by Conrad [Ref. 8] and 
ШЕГІ |Кет. 9| апа many others. The results have all 
indicated that the flow is pressure dependent as follows. 
In zone I the tube was collapsed and there was no flow. In 
zone II the flow was proportional to Р, - Pas and the tube 
was observed to neck near the downstream end. In zone III 
the flow was unobstructed and proportional to Р, - P5. 
These results were obtained using a latex rubber tube as a 


model surrounded by an airtight box. The fluid involved 


was water. The results were based on a single length-to- 
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diameter (L/D) ratıo and a single tube tensıon in each 
study. 

Ihis experiment was conducted using air аз the fluid 
Border To verify the previous results and to study the 
effects of variations in L/D ratio and tension in the tube 


wall. 


B. EXPERIMENTAL APPARATUS 

The experimental apparatus used in this investigation 
was designed and constructed similar to the previous models 
of Conrad [Ref. 8] and Holt [Ref. 9]. Alterations were 
miner Oduced where necessary to simplify the acus and 
cxt 3t to air flow. The system is shown in Figures 12 
апа 15. The details of the equipment were as follows. 

An airtight box was constructed with dimensions 14" by 
10" by 10" using 1/2" Plexiglas sheet. The inlet to the 
box was a decreasing radius nozzle constructed from powdered 
aluminum-filled epoxy resin. It converged smoothly from an 
outside diameter of 5-1/2'' to an inner diameter of 1" 
dekuminum pipe. The outside diameter of the pipe was 1-1/4". 
The exit from the box was via the same type pipe. The mouths 
Ене two pipes меге separated by 10", and radiused to 
ensure smooth flow. Variations in the L/D ratio were 
attained using several lengths of Plexiglas tubing adjusted 
Nereo ape lNnemmnstdesatanceen Of the 
Plexiglas БЫА was 1-1/4", and the outside diameter was 


1-1/2". The ШОТ of each tube was also radiused. A I" 
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gate valve was mounted between the exit and the vacuum 
source in order to adjust the pressure. Pressure in the box 
was.regulated by a 3/4" bleed from the exit pipe (taken 
upstream of the 1" gate valve) equipped with a 3/4" gate 
valve. Static pressure ports were located at the inlet, 
exit, and top of the box. A pitot tube was inserted in 

the exit flow. The pressure ports were connected to an 
inclined water manometer. The suction source was a Singer 
vacuum cleaner, Model D-5. The collapsible tube was a thin 
rubber prophylactic with a limp length and diameter of 7" 
and 1-3/8", respectively. It was stretched over the 
Plexiglas tube at each end, and fixed by rubber O-rings and 


bands. 


O EXPERIMENTAL PROCEDURE 
The basic experimental procedure consisted of variation 
a m tube tension for each of three L/D ratios.. 
Ihe 1" gate valve was opened to adjust the inlet pressure 
Р, to a value of 2.0 in H, 


maintained at this level for all investigations. One of 


O below atmospheric pressure, and 


three L/D ratios was then chosen and the tube tension was 
adjusted for each of three values. A new L/D ratio was 
chosen and the procedure was repeated. The pressure in the 
box Р was varied from 5.0 in Н, 


collapse of the tube. As a tube was stretched, its diameter 


O below atmospheric up to 


decreased, and the tension of the tube was quantified by 


the pressure difference Р, - Р, necessary to maintain 
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кога ет Ес та тнекпре [СиБе walls parallel to the pipe 


walls at the same diameter). 


D. RESULTS AND DISCUSSION 

The experimental study of air flow in a collapsible tube 
puuucedosignaficant results. The data collected are tabu- 
lated in Appendix B and graphically shown in Figures 14, 15, 
wao. In general, the results obtained were similar to 
those of previous work with water, and additional new infor- 
mation was observed. 

In zone I the flow was discovered to differ from Ze nor 
depending on the initial tension of the tube. In the low 
tension case, flow was found to exist for Pa greater than 
Р, Ше гот со аре әс пе tube. Generally, the effect of 
decreasing tension was to increase the pressure Pa attained 
Before collapse. In zone II the slope of the dynamic pres- 
mime curve (directly related to the mass flow) rapidly in- 
creased as Р, passed P^; and leveled as Р, approached (and 
17 In zone III the major deviation 
from previous results occurred. As the difference P, - P 


1 2 


was maintained as a constant in each case, the dynamic pres- 


passed in some cases) P 


sure should have remained constant for all Р, less than Р,. 


The results of this investigation indicated that the flow 
was constant until parallel flow was attained. However, 


flr ther decrease in P, caused the cross section to balloon, 


A 


accompanied by a drop in the dynamic pressure. 
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де depuccede s ii MI ааа c the effect of 
varying L/D ratios Was tomincwease the Р, possible before 
 ИТарве, vith decreasing L/D. 

Three qualitative phenomena were observed during 
œp r imentation. First, just prior to collapse, the tube 
walls developed longitudinal striation$ near the center of 
the tube. Second, collapse was a highly non-steady phenome- 
ПОП. It was not complete, steady collapse, but a rapid 
oscillatory motion of the tube walls from complete collapse 
emene Striated State described above. Finally, the col- 
lapse was observed to take place at a position approximately 
1/5 of the tube length from the downstream end for each 
value of L/D. The "ballooning" effect was also observed to 


begin at the downstream end. 


fee CONCLUSIONS 

The results of the investigation of air flow in a thin- 
ШЕШ ес collapsible tube justified the following conclusions: 

1. The experimental apparatus used was an accurate 
model of flow in a collapsible tube utilizing new parameter 
variations. 

2. The length-to-diameter ratio of the flexible tube 
was found to have a significant effect on pressure-flow 
relationships. 

5. Ibe longitudinal tension of the flexible tube was 
also found = have те ето сапе естеств on the 


pressure flowmreladtionships. 
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FIGURE 3.  ONE-DIMENSIONAL MODEL OF THE LUNG 
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BIGURE 15. EXPERIMENTAL APPARATUS 
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FIGURE 15. FLOW VARIATION WITH ALVEOLAR PRESSURE, 40 = 4.0. 
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FIGURE 20.  UNSTEADY COLLAPSE. 
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APPENDIX A 


NUMERICAL SOLUTION DATA 


N, eonegenrrasıon. durıne Inspiration, T.. = 


in 
Distance Concentration (с/с) 

(х/%) N = 6 N = 8 i N = 10 
Zi PRU .000 .000 .000 
-4.0 .000 .000 ‚000 
E .006 . 002 .000 
2) ‚050 ‚050 ОТЕ 
= МЕА SO „103 ‚087 
-1.0 . 216 .193 | E 
-0.4 . 389 2567 SO 

0.0 ‚545 ‚549 .554 

0.2 1.000 1.000 1.000 

0.6 1.000 1.000 1.000 

0.98 1.000 1.000 1.000 
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N, Concentration during Inspiration, N= 10 


ies tance Concentration (с/с) 

(x/2) цэ. Di Me T,, * 48 E. 
-5.0 . 000 .000 | . 000 .000 
-4.0 ‚000 ‚000 „000 „000 
-3.0 ‚000 ‚000 ‚000 ‚000 
-2.0 0.17 ‚000 . 000 .000 
-1.4 ‚081 ‚000 ‚000 „000 
-1.0 E? . 000 . 000 . 000 
BNET .387 „000 | „000 „000 

0.0 „554 ‚000 222-2000 „000 

(1 .999 2 . 870 . 301 21128 

0.2 1.000 „998 ‚641 „521 

0.4 1.000 1.000 .970 „746 

0.6 1.000 1.000 .999 .966 

0.8 1.000 1.000 1.000 .999 


uS 1.000 1.000 1.000 1.000 
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N, Concentration during Expiration, N =710 


Distance Concentration (c/co) 
(х/%) Т. =12 17-28 90 
lex 14 lex 12 To Төх 12 тте 2225 
-5.0 2592 ЕП) . 586 20723 .344 .642 
-4.0 559 185 595 080 ST 646 
ЕО -980 . 201 ‚601 ‚087 2 ‚646 
xn) ‚610 2216 ‚609 .095 05 ‚656 
-1.0 20:57 m Sl ‚616 „102 eS ae. „661 
0.0 2052 ‚ 246 2029 „110 „380 . 666 
0.1 205 .439 . 7 20 ‚216 „479 .728 
022 ‚984 ЕЭ nonu .418 ‚650 -929 
03 .999 .898 „945 ‚628 29D .904 
0.4 1.000 2375 „982 . 800 ‚894 „949 
0.6 1.000 „999 998 .969 217: .988 
0.8 1.000 127000 1.000 „999 „998 2996 
0.98 1.000 1.000 1.000 1.000 1.000 „ООО 
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ШЕ а. 





Washout Data, T. = 12 
in 


Tas Concentration (с/с) 
(ugt/£) М = 6 N = 12 
2.4 ‚000 | ‚000 
4.8 SOT .061 
EY SN 298 
9.6 “НЭЭ . 438 
72.0 . 540 „550 
14.4 ‚603 е 
16.8 2051 | ‚650 
E ‚689 2092 
2170 20221 200 
24.0 .748 ВУ БІВ 


SU 





12.0 
18.0 
24.0 
S0. 0 
50.0 
42.0 
48.0 
54.0 
50.0 
66.0 
2.0 


Washout Data, 


ЫМ - 6 


3% 


D С 
in 


Conceneration- (c/c 


0: 


N 


= IZ 





ЦЕ? 0 
24.0 


36.0 


Washout Data, T._ = 96 


N 


zo 


ЭЭ 


11 


Concentration (c/c 


07 


п 


Н 
— 
N 


2507 
MOTA 
‚687 
НИВИ 





APPENDIX E 


THIN-WALLED COLLAPSIBLE TUBE DATA 


All pressures are inches of water below atmospheric 


Pressure. Р, СЕЛ |, огт э!? cases, 


 ШШ- 55 АР = 2.3 


р 
2 Pa q 
2.85 5.0 Е 
211 4.5 ‚60 
2.85 4.3 760 
21 4.0 55 
2.88 2 | ‚50 
2.85 570 ‚45 
2.85 2.8 ‚45 
2.85 2.6 ‚40 


L/D - 5.5, AP, - 1.9 


2.80 3.0 2605 
ZO e 4.5 . 60 
2221) | 4.0 219,9 
2.80 Ээ? ОО 
2.80 22:15 229 
2121! 221! . 50 
221 2.6 ‚45 
2.80 2.4 ‚40 
2.80 26123 . 40 
22 800. 22 22.35 
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Э 
Э 
5 
го 
2225 
nu 
215 
Э 


2 ЕЭ 


‚60 
‚860 
. 80 
. 80 
. 90 
‚60 


L/D 


L/D 


sr 


4. 


ЫП 
р 


55 


$05 
.65 
ОЮ 
ае 
255 
-50 
. 45 
. 40 
.40 


.60 
255 
250 
‚45 
„45 
.45 





. 80 
.80 
. 80 
. 80 
.90 
„80 
‚860 
„80 
. 80 


ЕЭ 
SiS 
29 
D 
22125 

До 
22 
БЭЭ 
TUS 
sho 
2025 
el 


L/D 


4. 


‚60 
. 50 
.45 
.45 
„45 
‚40 
о 
155 


295 


0 
„45 
„40 
„40 
„40 
„40 
„40 
.40 
„35 
‚50 
ЕТ 
‚30 


ша 








2.5 
5 
15 
ES 
„5 
То 
S 
m5 
SES 


210 
-70 
ЕГО 
270 
411! 
20 
D 
‚70 
‚70 
400 
жи 


L/D 


zs 


37 


2450 
‚45 
‚40 
. 40 
. 40 
.40 
. 40 
ЖЭЭ 


229 


250 
. 45 
‚40 
‚40 
‚40 
‚40 
‚40 
‚40 
235 
229 


m 





255 
. 65 
40:5 
205 
205 
#05 
305 
-09 
‚ 65 
ӨБ 
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‚50 
‚ 45 
. 40 
. 40 
99 
223 
FI 
‚50 
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